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ABSTRACT 

We present an optical nebular spectrum of the nearby Type la supernova 201 Ife, 
obtained 981 days after explosion. SN2011fe exhibits little evolution since the +593 
day optical spectrum, but there are several curious aspects in this new extremely 
late-time regime. We suggest that the persistence of the ^ 5800 A feature is due to 
Nal D, and that a new emission feature at ~ 7300 A may be [Call]. Also, we discuss 
whether the new emission feature at ^ 6400 A might be [Fel] or the high-velocity 
hydrogen predicted by Mazzali et al. The nebular feature at 5200 A exhibits a linear 
velocity evolution of ~ 350 km s“^ per 100 days from at least +220 to +980 days, but 
the line’s shape also changes in this time, suggesting that line blending contributes 
to the evolution. At ^ 1000 days after explosion, flux from the SN has declined to a 
point where contribution from a luminous secondary could be detected. In this work 
we make the first observational tests for a post-impact remnant star and constrain its 
temperature and luminosity to T > 10^ K and L < 10'* Lq. Additionally, we do not see 
any evidence for narrow Ha emission in our spectrum. We conclude that observations 
continue to strongly exclude many single-degenerate scenarios for SN 201 Ife. 

Key words: supernovae: general — supernovae: individual (SN 201 Ife) 


1 INTRODUCTION 

Supernovae of Type la (SNela) are thermonuclear explo¬ 
sions of carbon-oxygen white-dwarf stars. As standardisable 
candles they are a very valuable tool for modern cosmology, 
but their progenitor scenario(s) and explosion mechanism(s) 
are not yet well understood (e.g., Howell 2011). For the for¬ 
mer, the two leading candidates are the double-degenerate 
scenario, in which the primary merges with another carbon- 
oxygen white dwarf or accretes matter from a helium white 
dwarf, and the single-degenerate scenario, in which matter 
is accreted from a red-giant or main-sequence companion. 
Most of the current efforts toward identifying the progeni¬ 
tor are focused on pre- or post-explosion imaging to directly 
detect the progenitor star and/or its companion (e.g., Li et 
al. 2011; Kelly et al. 2014), and looking for evidence of the 
SNIa ejecta interacting with circumstellar material (CSM) 

* E-mail: melissagrahamQberkeley.edu 


— the amount, composition, and distribution of which the¬ 
oretically depends on the progenitor system (e.g., Leonard 
2007; Patat et al. 2007; Sternberg et al. 2014). 

SN 201 Ife was discovered by the Palomar Transient Fac¬ 
tory (Ran et al. 2009; Law et al. 2009) within hours of its 
explosion on 23.7 Aug. 2011 (dates are UT throughout this 
paper) and spectroscopically classified as a normal SNIa 
(Nugent et al. 2011). The host galaxy is MlOl (NGC 5457), 
just 6.4 + 0.7 Mpc away (Shappee & Stanek 2011) with a re¬ 
cession velocity of Vrec = 241 km s“* (de Vaucouleurs et al. 
1991). Richmond & Smith (2012) present photometric mon¬ 
itoring, reporting that SN 2011fe reached peak brightness on 
11.5 Sep. 2011. Nugent et al. (2011) used the very early-time 
light curve to constrain the size of the progenitor to white- 
dwarf scale (later refined by Bloom et al. 2012), and they 
show that SN 201 Ife exhibited no evidence of interaction be¬ 
tween the SN ejecta and the companion star as predicted by 
Kasen (2010). Li et al. (2011) use pre-explosion Hubble Space 
Telescope (HST) archival images to exclude bright red giants 
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and most main-sequence helium stars as mass-donor com¬ 
panions to the white dwarf. Graur et al. (2014) show that 
archival narrow-band HST images centred on Hell A4686 
indicate no supersoft X-ray source existed at this location, 
contrary to expectations for the single-degenerate scenario. 
These studies eliminate many single-degenerate companions 
for SN2011fe. 

The proximity of SN 201 Ife allows it to be monitored to 
extremely late times (e.g., Kerzendorf et al. 2014; Tauben- 
berger et al. 2014), and in this evolved state two tests of 
the progenitor system become feasible. First, the lower lu¬ 
minosity of the SN la may allow for even a relatively small 
mass of hydrogen in the system to be observed via Ha emis¬ 
sion (e.g.. Marietta et al. 2000; Leonard & Filippenko 2001; 
Mattila et al. 2005; Leonard 2007; Lundqvist et al. 2013; 
Shappee et al. 2013a; Lundqvist et al. 2015). Hydrogen is 
expected to be present if the companion was a red giant or 
main-sequence star, but not in any appreciable amount if 
the companion was a carbon-oxygen or helium white dwarf. 
Second, the lower luminosity of the SN la may allow for a di¬ 
rect detection of the remaining companion star in the single¬ 
degenerate scenario, because the secondary is theoretically 
predicted to increase in temperature and luminosity in the 
years after impact with the SNIa ejecta (e.g.. Pan et al. 
2013; Liu et al. 2013; Shappee et al. 2013b). 

In this work, we use an optical nebular spectrum ob¬ 
tained 981 days after explosion to constrain the amount of 
hydrogen present in the system, and the temperature and lu¬ 
minosity of a putative post-impact remnant companion star. 
In ij^we present our observations, and in i|^we analyse our 
spectrum, comparing it to previous nebular observations in 
2013 and 2012. Section [3.11 demonstrates a nondetection of 
Ha and places an upper limit on the mass of hydrogen in 
the system, and §3.2| provides two independent tests for the 
signature of a post-impact remnant star (PIRS) in our spec¬ 
trum of SN 201 Ife. We summarise and conclude in 


2 OBSERVATIONS 

We obtained spectra of SN 2011fe with the Low Resolution 
Imaging Spectrometer (LRIS; Oke et al. 1995) at the Keck 
Observatory on 2014 May 01 — that is, 981 days after explo¬ 
sion (we use observer-frame time). At such a late phase, the 
source was too faint to appear in the LRIS guide-camera im¬ 
ages recorded during spectroscopy, so ~ 3.5 weeks later we 
obtained images in the g and Rs bands with LRIS on 2014 
May 25 (-1-1006 days). Here we present and discuss first the 
imaging in § |2.1| and then the spectroscopy in § |2.2| All of 
the spectra are available on WISeRElQ (Yaron & Gal-Yam 
2012 ). 

2.1 Imaging 

With Keck-|-LRIS, we acquired three images in the g fil¬ 
ter of 180, 220, and 230 s duration (increasing the expo¬ 
sure time because the red-side readout time is significantly 
longer than that of the blue side), and three images in the 
Rs filter of 180 s each. SN2011fe is in a crowded field with 

^ http://wiserep.weizmann.ac.il 



Figure 1. Images from Keck+LRIS on 2014 May 25 in the LRIS 
g and Rs filters (left and right, respectively). Pixels have been 
grey-scaled to show the faint SN 201 Ife (circle) in a field crowded 
by brighter sources and contaminated by diffuse background light 
from MlOl. The blue box represents the slit width and orienta¬ 
tion used for the Keck+LRIS spectrum obtained on 2014 May 
1, and the size of the circles is set equal to the full width of the 
aperture used to extract the 1-dimensional spectrum. Despite the 
crowded field, our slit width, orientation, and chosen aperture 
limit contamination from nearby stars. 

much host-galaxy background light, and Sloan Digital Sky 
Survey (SDSS) standard stars are not available within our 
small field of view. For photometric calibration, we also ob¬ 
tained two 1 s images in each of two regions in the GFHT 
Legacy Survey Deep Field 3 immediately after SN 2011fe; 
these fields were chosen because they are well covered (e.g., 
by SDSS) and are only ~ 2° from SN2011fe. We applied 
typical image reduction and analysis procedures to all fields, 
and used SOURCE Extractor (Bertin & Arnouts 1996) to 
find and measure source photometry (using aggressive de¬ 
blending parameters for the crowded field of SN 201 Ife). Our 
final stacked images are shown in Figure 

On 2014 May 25, we had an LRIS blue-side shutter 
failure and the chip remained illuminated during readout. 
The blue-side readout time is 42 s, so this affects the Deep 
Field images used to determine our zeropoint more than it 
affects the SN20life images. We use SDSS star magnitudes 
from Data Release 10 (Ahn et al. 2014) to calculate a sim¬ 
ple zeropoint {zg — mg^sDSS — nig^ram) and find that this 
increases by ~ 0.2 mag across the chip. We take this into 
account by using the mean zeropoint, which also represents 
the zeropoint at the centre of the chip where SN 20life is, 
and adding 0.2 in quadrature to the uncertainty. We esti¬ 
mate that the brightness at the position of SN 2011fe was 
rrig ~ 24.3 ± 0.3 mag on 2014 May 25. The line-of-sight 
Galactic extinction is much smaller than our uncertainty, 
Av = 0.023 mag from the NASA Extragalactic Database 
(NEEQ Schlafiy & Finkbeiner 2011; Schlegel, Finkbeiner & 
Davis 1998). For comparison, Kerzendorf et al. (2014) find 
that rUg = 23.43 ± 0.28 mag on 2014 March 7; assuming a 
decline rate of 0.01 mag per day, the predicted brightness 
for our observation is rUg « 24.0 mag. Given our large un¬ 
certainties, our results are not discrepant with this predic¬ 
tion — however, there is evidence that the rate of decline 
is probably slower than 0.01 mag per day at this very late 
phase (e.g., Kerzendorf et al. 2014). On the red side, we es- 

^ http://ned.ipac.caltech.edu/ 
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timate mu « 24.6 mag with similarly large uncertainties. As 
we do not use this in our analysis, we leave a more accurate 
photometric measurement for a future publication, includ¬ 
ing template image subtraction to remove background light. 
Note that Kerzendorf et al. (2014) also have not subtracted 
a template image. 

In §3.2[ our constraints on the presence of a PIRS de¬ 
pend on an accurate flux calibration for our spectrum, es¬ 
pecially at the blue end, and our results are more conser¬ 
vative if we adopt a brighter y-band magnitude. Synthe¬ 
sised spectrophotometry from our reduced Keck-|-LRIS spec¬ 
trum (described below) indicates mg ~ 23.9 ± 0.2 mag and 
Rs > 24.2±0.5 mag, in agreement with (but slightly brighter 
than) our photometry. Since this is quite sufficient for the 
analysis, we make no further adjustments to the flux cali¬ 
bration of the spectrum. 

2.2 Spectroscopy 

For our LRIS spectra we used the 1.0” slit rotated to the 
parallactic angle (Filippenko 1982), which was 200° east 
of north, to minimise the effects of atmospheric dispersion 
(in addition, LRIS has an atmospheric dispersion correc¬ 
tor). The size and orientation of the slit is shown as the 
blue box in Figure An offset star was used for target ac¬ 
quisition because SN2011fe was too faint to appear in the 
guider camera. In our LRIS configuration, coverage in the 
blue with the 600/4000 grism extends over A = 3010-5600 A, 
with a dispersion of 0.63 A pixel”^ and a full width at half¬ 
maximum intensity (FWHM) resolution of ~ 4 A. Coverage 
in the red with the 400/8500 grating extends over A = 4500- 
9000 A, with a dispersion of 1.16 A pixel”^ and a resolution 
of FWHM « 7 A. We did not dither along the slit between 
exposures, instead staying at the same CCD location in all 
four 1200 s exposures. 

These spectra were reduced using routines written 
specifically for LRIS in the Carnegie Python (CarPy) 
package. The two-dimensional (2D) images were flat fielded, 
corrected for distortion along the y (slit) axis (i.e., recti¬ 
fied), wavelength calibrated, cleaned of cosmic rays (imper¬ 
fect), and stacked with a median combine before defining the 
apertures and extracting the one-dimensional (ID) spectra 
(sky subtraction is discussed below). The width of the ex¬ 
traction aperture was ±5 pixels (±0.68”) on the blue side 
and ±4 pixels (±0.54”) on the red side, as shown by the 
radius of the green and red circles in Figure The seeing 
is estimated to have been ~ 0.9”(FWHM) during the obser¬ 
vations. Although the field is crowded. Figure]^ shows that 
the nearest neighbours do not contaminate our spectrum. 
Blue and red standard stars were used to determine the sen¬ 
sitivity functions and flux calibrate the data for each side, 
respectively. 

The input parameter space for these routines was care¬ 
fully explored to ensure the best possible results, and we 
found that the method for removing the sky background had 
the largest influence on the final spectrum. The two options 
for treating the sky background are as follows. “Redux A”: 
create a sky frame for each 2D image and subtract it prior 
to median combining the images and performing aperture 
extraction (Kelson 2003). “Redux B”: define background 
apertures on either side of the trace during aperture ex¬ 
traction from the stacked 2D image and subtract from the 
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Figure 2. The reduced and calibrated LRIS spectra of SN 2011fe, 
showing the small systematic difference introduced by the treat¬ 
ment of the sky background for this faint object. For “Redux A” 
(blue) we created a sky frame for each 2D image and subtracted it 
prior to the median combine and aperture extraction; for “Redux 
B” (orange) we defined background apertures on either side of 
the trace during aperture extraction from the median-combined 
2D image. We show the resulting spectra for the whole wave¬ 
length range (top) and the region around Ha only (bottom), in 
which the solid red line marks the expected position of Ha at the 
recession velocity of SN 2011fe. Spectra are binned as described 
on the plot; horizontal lines at zero flux are included to guide 
the eye. The excess noise at 5600 A is at the dichroic’s crossover 
wavelength. 


ID spectrum (i.e., set the background parameter to “fit” 
in IRAF task apall). We used 4-pixel apertures adjacent to 
the extraction aperture. The results of these two reduction 
versions are compared in Figure 

Sky subtraction is the largest (although still small) 
source of systematic uncertainty. On the blue side (A < 
5500 A), the continuum of “Redux A” does not go to zero 
flux at A « 4900 A like the continuum of “Redux B.” This 
is not necessary, as both residual photospheric emission and 
low-level flux from emission lines can produce some contin¬ 
uum flux which is actually from the SN. However, it might 
indicate that an incomplete sky subtraction is achieved on 
the blue side when a sky frame for each 2D image is created 
and subtracted. In Figure[^we can see that the choice of sky 
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subtraction leads to a small difference in relative emission¬ 
line strengths on the blue side, but this will not affect our 
particular analysis. On the red side, the sky background 
is concentrated more into emission lines, and has less flux 
in the continuum, than on the blue side. On the red side, 
however, “Redux A” has a lower scatter and fewer residual 
spikes from the sky emission lines because the 2D subtrac¬ 
tion is more sophisticated than the ID background aper¬ 
tures used in “Redux B.” We proceed using both versions 
of our reductions so that our conclusions are robust to sys- 
tematics introduced by the data reduction for this very faint 
SN. Where necessary, the following analysis will state which 
reduction version is being used. Before any analysis is per¬ 
formed, the spectra are deredshifted into the rest frame of 
the host galaxy, MIDI, and corrected for line-of-sight Galac¬ 
tic extinction E{B — V) ~ 0.008 mag (Schlegel, Finkbeiner, 
& Davis 1998). 


2.3 Earlier Nebular Spectra of SN 2011fe 

For comparison, we also use nebular spectra of SN 2011fe ob¬ 
tained under the Berkeley Supernova la Program (BSNIP; 
Silverman et al. 2012) in 2012 and 2013, shown in Figure 
The 2012 spectrum was taken on 23 Aug. 2012, 365 days 
after explosion, with the Kast spectrograph at Lick Obser¬ 
vatory (Miller & Stone 1993), using the 600/4310 grism and 
the 300/7500 grating, a long slit 2.0” wide, and a 1800 s ex¬ 
posure time. This spectrum is previously unpublished and 
will also appear in Mazzali et al. (2015). The 2013 spec¬ 
trum was obtained on 8 April 2013, 593 days after explosion, 
with DEIMOS (Faber et al. 2003) on Keck II, using the 600 
lines mm“^ grating, the GG455 order-blocking Alter, the 
1.0” slit of the long variable multislit (LVM) slitmask, and 
a 1200 s exposure time. This spectrum is previously unpub¬ 
lished. The data were reduced using common reduction and 
analysis procedures (e.g., Silverman et al. 2012). 

In Figurej^we have flux scaled the 2012 and 2013 spec¬ 
tra to the 2014 spectrum by multiplying the 2012 and 2013 
spectral flux by constant normalisation factors. The normal¬ 
isation factors are set to values that minimise the mean ab¬ 
solute flux residuals in particular wavelength regions (0.005 
and 0.122 for the 2012 and 2013 spectra respectively). These 
regions, shown in Figure]^ are chosen to be areas on the red 
side in which the feature shapes were well matched between 
epochs and minimal evolution is expected. The red side was 
also chosen for this minimisation because most of the anal¬ 
ysis will be done on the blue side. 

In pTT] we want to compare the variance of the 
2014 spectrum with that from 2013, but unfortunately the 
DEIMOS spectrum from 2013 does not extend bluer than 
~ 4500 A. As a compromise, since the region blueward of 
4500 A exhibits similar features in 2012 and 2014 (Figure]^, 
we supplement the 2013 spectrum with the 2012 spectrum in 
the range 3500-4500 A to create a “composite” 2012-1-2013 
SN 2011fe spectrum. First, the 2012 spectrum is shifted in 
wavelength such that the broad peaks at 4400 and 5200 A 
match with their positions in 2014 (~ 60 A). Then, the 2012 
and 2013 spectra are flux scaled to the 2014 spectrum as de¬ 
scribed above, and then joined to produce a final composite 
spectrum. The application of this composite spectrum is de¬ 
scribed in j ]3.2.1| 



Figure 3. A comparison of our spectrum of SN 2011fe from 2014 
(black) to nebular spectra from 2012 (pink) and 2013 (blue). The 
latter two have been flux scaled to the 2014 spectrum by deter¬ 
mining the normalisation factor that minimises the absolute value 
of their flux residuals in the wavelength ranges specified on the 
plot (0.005 and 0.122 for the 2012 and 2013 spectra respectively, 
as described in the text of § |2.3| no wavelength shifts are ap¬ 
plied here). The ordinate values therefore apply only to the 2014 
spectrum. The ionised elements commonly attributed to each line 
are labeled with coloured bars, and our proposed line identifica¬ 
tions are labeled in grey. While the new 2014 spectrum resembles 
the 2013 spectrum, both are significantly different from the 2012 
spectrum, particularly in the lines of [Felll]. 

3 ANALYSIS 

Early-time photospheric-phase optical spectra of SN 2011fe 
were analysed by Parrent et al. (2012), who classify it as 
spectroscopically normal and a member of the “low velocity 
gradient” group, referring to the change in velocity of the 
Sill A6355 line (Benetti et al. 2005). Mazzali et al. (2014) 
apply abundance tomography models to near-ultraviolet and 
optical spectra of SN 2011fe around maximum light and 
show that ®®Ni is distributed in an extended region; that 
SN 201 Ife probably had a lower density and larger volume of 
nucleosynthetic products; and that the amount of iron in the 
outer layers of SN 2011fe is consistent with subsolar metal- 
licity, Ziife « 0.5 Zq. A subsolar metallicity is also sup¬ 
ported by the comparative analysis of photospheric through 
nebular-phase spectra of SN 2011fe and SN 2011by by Gra¬ 
ham et al. (2015). The extensive collection of analyses of 
SN 201 Ife agree that it was a typical example of a SN la. 

In Figure we compare our 2014 spectrum with flux- 
scaled spectra of SN2011fe from 2013 and 2012. We find 
that over the last year there has been little evolution in 
the spectral features; the 2013 and 2014 spectra are very 
similar to each other, and different from the 2012 nebular 
spectrum. The similarity of the spectral features as the total 
flux declines is evidence of no (or negligible) contributions 
from a light echo, consistent with previous findings that 
SN 2011fe exploded in a “clean” environment (e.g., Chomiuk 
et al. 2012; Horesh et al. 2012; Patat et al. 2013; Pereira et 
al. 2013). We point out the following changes between the 
2012 and 2014 spectra: (1) the decline of the [Felll] line at 
A4700 A, (2) the feature at ~ 5800 A where, at earlier times, 
[Coin] appears, (3) the potential rise of a new broad fea- 
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ture at 6400 A, just blueward of where the [Co III] ~ 6600 A 
feature declined after 2012, (4) a decline in the [Fell] and 
[Nill] feature at ~ 7300 A, and (5) a general redshifting 
of all spectral features. These changes are discussed below 
(and in some cases also by Taubenberger et al. 2014). 

The decline of doubly ionised species [Fe III] and [Co III] 
can be attributed to the decreasing temperature in the rem¬ 
nant ejecta material. Curiously, the feature at 5800 A, com¬ 
monly thought to be [Co III] , appears quite strongly in both 
our 2013 and 2014 spectra, but is not seen by Taubenberger 
et al. (2014). We think this discrepancy is caused by a co¬ 
incidence of this feature with their spectrograph’s dichroic 
at 5750 A, as we cannot attribute this feature to be an ar¬ 
tifact of our reduction pipeline (it appears quite plainly in 
our raw, sensitivity- and flux-calibrated data, for all versions 
of our reduction process). Given the decline of the [Co III] 
6600 A feature, another species is probably responsible. [Fe I] 
is a likely culprit, and it has plenty of transitions in this re¬ 
gion, such as 5640, 5695, and 5709 A (Smith & Wiese 1973). 
However, this line could also be Nal D, which has long been 
predicted at this location by models, and first considered 
(but rejected) in nebular spectra by Kuchner et al. (1994). 
Dessart et al. (2014) present a review of the literature for 
and against Nal D at > 100 days (see their Section 2), and 
show that [Coin] is the main contributor for the first two 
months after peak brightness, but at later times the Na may 
finally be visible because it does not decay like ®®Co. For ex¬ 
ample, Nal D has been attributed to the 5800 A feature for 
the normal SNIa 2003du by Tanaka et al. (2011). 

Regarding the rise of the broad feature at 6400 A, 
Taubenberger et al. (2014) attribute this to [Fel] or [Ol], 
pointing out that detailed spectral modeling is required to 
confirm it. We discuss in § |3.1.1| the possibility that this is 
broad, high-velocity emission from hydrogen. Further red- 
ward, we see that the neighbouring [Fell] A7155 and [Nill] 
A7370 features have declined, and that this region is now 
dominated by a third peak at ~ 7300 A, which could be the 
redshifted line of [Fell], or [Call] AA7292, 7324. Emission 
from Na and Ca may suggest the presence of circumstellar 
material from the progenitor system. 

We find that the continuous shift in the nebular fea¬ 
tures toward redder wavelengths is part of a linear trend 
that starts by at least 200 days after peak brightness. We 
measure the central wavelength of the [Fe II] 5200 A feature 
in the 2014, 2013, and 2012 nebular spectra used in this 
work, as well as in the -1-300 and -1-200 day spectra from 
Graham et al. (2015), by smoothing the spectra and find¬ 
ing the wavelength at the centre of the FWHM. We find a 
linear evolution of ~ 6 A (~ 350 km s“^) per 100 days for 
SN 201 Ife. This is somewhat lower than the typical evolution 
found by Silverman et al. (2013) for a larger sample of SNIa 
nebular spectre]^ for which the average velocity of the neb¬ 
ular features decreases by 1000-2000 km s“^ per 100 days. 
They also find that the velocity appears to evolve to u « 0 
km s“^ by ~ 300 days past maximum brightness for the 
SNela in their sample (e.g., see their Figure 4), suggesting 


^ Note that Silverman et al. (2013) use the velocity of the [Felll] 
~ 4700 A feature, whereas we must use the [Fell] ~ 5200 A 
feature, because the former does not appear in our 2013 or 2014 
spectra. 


that continued velocity evolution until -flOOO days results 
in a redshift. Taubenberger et al. (2014) also interpret their 
2014 spectrum of SN2011fe as exhibiting a redshift in the 
rest frame. However, for SN2011fe, the 5200 A feature is 
not symmetric, and we see that the wavelength of the peak 
flux moves from the blue to the red side of the feature over 
time. Instead of using the wavelength at the centre of the 
FWHM, if we use the wavelength at peak flux to calculate 
the linear evolution in the line velocity, we find ~ 11 A 
(~ 640 km s“^) per 100 days for SN2011fe (i.e., the line’s 
peak moves more rapidly redward than the line’s centre). 
The fact that the line morphology changes significantly sug¬ 
gests that line blending cannot be ruled out as the source of 
the apparent shift, and we conclude that we cannot be sure 
that the nebular features are exhibiting a redshift in the 
rest frame. We hope that future spectral modeling and/or 
observations of extremely late-time SNela will provide ex¬ 
planations for all of these observations. 

In 2014, the flux of SN 201 Ife has declined far enough 
that tests of the progenitor companion using a nebular-phase 
spectrum are possible. In the following sections we perform 
three tests to constrain the properties of a putative compan¬ 
ion star. 

3.1 An Upper Limit on the Mass of Hydrogen in 
the Progenitor System of SN 2011fe 

Type I supernovae are classified by a lack of hydrogen in 
their spectrum, but some progenitor scenarios, such as mass 
accretion from a main-sequence or red-giant star, predict 
that some hydrogen can be seen. For example, hydrogen 
could be stripped from a main-sequence star and be moving 
with < 1000 km s“^ (dA « 22 A; Liu et al. 2013, Figure 9). 
Nugent et al. (2011) found that there was no main-sequence 
or red-giant companion along our line of sight to SN 2011fe, 
as the interaction between SN ejecta and such a star would 
increase the blue flux at very early times (Kasen 2010), but 
a main-sequence star might have existed on the far side. 
If hydrogen is present in the system, we might be able to 
see it at late times - and if we do, then hydrogen cannot 
be ruled out as being absent in any SN la that has only 
early-time nondetections of hydrogen. We have carefully and 
thoroughly scrutinised our 2D reduced spectra but find no 
trace of an emission line at or anywhere near the position 
of SN 2011fe, before or after background subtraction. (We 
do see host-galaxy emission in the lower half of the LRIS 
detector, on the other side of the chip gap from SN 2011fe.) 

As shown in the bottom plot of Figure [^ there is no 
obvious emission seen in the ID spectrum either. To quantify 
the maximum potential equivalent width of an Ha emission 
line “hiding” in our data, we follow a similar analysis to 
that of Shappee et al. (2013a), who constrain the amount of 
hydrogen in the system of SN 20life using their deep spectra 
taken up to 275 days past maximum light. 

Marietta et al. (2000) model the evolution of hydro¬ 
gen from the companion star to 380 days after maximum 
brightness. Applying this model to SNIa 200lei (which 
was normal, like SN2011fe), Mattila et al. (2005) find that 
M = 0.05 Mq of material would create an Ha line with 
peak luminosity ~ 3.36 x 10^® erg s A“^. Shappee et al. 
(2013a) find that for the distance and Galactic reddening of 
SN 2011fe, this would be an observed peak of 6.71 x 10“^^ 
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Figure 4. The ID spectrum from our reduction version “A” (see 
and Figure]^, shown in the rest frame of SN2011fe in the 
region of Ho. Our LRIS spectrum from 2014 is given in black, 
the continuum from polynomial smooth in blue, and a Gaussian 
profile representing the maximum possible Hq feature in red. The 
Gaussian has a width of 22 A and height equal to 3 times the la 
RMS fluctuation around the continuum. (Note that the suspicious 
narrow line has a FWHM of 1.8-3.1 A, < 200 km s“^; it is likely 
to be noise.) 

erg cm“^ A“^. An expression for the maximum equiv¬ 
alent width of an emission line “hiding” in a spectrum, 
VFA(3cr), is derived and presented by Leonard Sz Filippenko 
(2001) and Leonard (2007): 

Wx{3a) = 3 AlVWline AX, (1) 

where A7 is the la root-mean square (RMS) of the fluctua¬ 
tions around the normalised continuum, Wiine = 22 A, and 
AX = 0.61 A (the resolution or bin size). To measure A7, 
we apply a least-squares polynomial smoothing filter (Sav- 
itzky & Golay 1964) with a width of 163 pixels (100 A), 
subtract the data from the smoothed continuum, and find 
the square root of the mean of the squared residuals in the 
range 6480-6700 A. The value of A7 is then this RMS fluc¬ 
tuation divided by the continuum. 

In !|^ we discuss our two reduction versions, “A” and 
“B,” and Figure shows that “Redux A” has less noise 
in the region around Ha owing to a more sophisticated 
background subtraction. With “Redux A,” we find that 
the maximum equivalent width of hydrogen allowed by the 
data is WA,A(3cr) ~ 24.0 A (for reduction version “B” the 
result is ITa,s(3ct) « 42.9 A, so our uncertainty here is 
AITa.a « 18 A). This is shown in Figure]^ 

An emission line of peak 7 = 6.71 x 10“^^ 
erg s“^ cm“^ A“^ and FWHM = 22 A on top of our 
data would have a whopping equivalent width W\{Mu = 
0.05 M 0 ) « 4000 A. If we scale linearly between the equiv¬ 
alent width and mass of hydrogen (as do Shappee et al. 
2013a), then our upper limit on the amount of hydrogen is 

Mh < X 0.05 Mo « 3 X 10"'‘ Mo, (2) 

4000 A 

with an uncertainty of AMh ~ 2.4 x 10“^ Mq. Our results 


are very similar to the limit produced by Shappee et al. 
(2013a), A7h = 5.8 x 10“"^ Mq, from a longer integration at 
an earlier phase. This uncertainty includes measurement er¬ 
rors only, and not systematics from model assumptions (see, 
e.g., Lundqvist et al. 2015). At this time, our estimate of Mh 
in Equation]^ is essentially a lower limit on the upper limit 
of hydrogen in the system, and our main result is simply the 
qualitative nondetection of Ha. 

3.1.1 Broad Ha at 64 OO A? 

Taubenberger et al. (2014) attribute the new broad feature 
at 6400 A (seen in Figures 1^ and 1^ to [Fel] or [Ol], and 
mention that detailed spectral modeling is required to con¬ 
firm this. We point out that Mazzali et al. (2014) argue that 
the presence and strength of carbon at early times indicates 
there may be some hydrogen (~ 0.01 Mo) in the outer lay¬ 
ers of SN 2011fe, which might be seen at high velocity (for 
SNela, “high velocity” refers to 15,000-20,000 km s“^). Hy¬ 
drogen at ~ 15, 000 km s“^ would appear at ~ 6250 A, and 
given its location on the surface of the SN la — not out 
in the CSM as considered in the analysis above — may be 
seen as a broad emission line. If the broad feature we see 
at ~ 6400 A is hydrogen, its velocity is ~ 7500 km s 
significantly lower than predicted. To estimate the equiva¬ 
lent width of the 6400 A feature in our 2014 spectrum, we 
aggressively bin (51 pixels, 32 A) and numerically integrate 
to find EW ~ 700 A. By Equation]^ this implies a mass of 
hydrogen of Mh ~ 0.009 Mo. Although this highly uncer¬ 
tain estimate is in good agreement with the predictions of 
Mazzali et al. (2014), the caveats described above regarding 
the conversion from equivalent width to mass apply here 
as well. Encouragingly, redshifted emission from hydrogen 
on the far side of the explosion may be seen at ~ 6700 A 
in Figure but upon further inspection of the 2D images 
we find a diffuse extended artifact at that location, so this 
interpretation is not secure. 

3.2 Limits on the Luminosity and Temperature of 
a Post-Impact Remnant Star Companion 

In the single-degenerate scenario, when the companion star 
of an exploding white dwarf has its outer layers blown away 
by the supernova, the PIRS becomes hotter and more lu¬ 
minous. Models for nondegenerate He-burning companions 
(Pan et al. 2012, 2013) have shown they reach thousands of 
Lo after 1-3 yr, and tens of thousands of Lo after 10 yr. 
These PIRS can be modeled as blackbodies (BBs), and 
given their high temperatures they are significantly blue. 
At ~ 1 yr after peak brightness, nebular emission from 
SN 2011fe would swamp any contribution from a companion, 
but by ~ 3 yr the contribution of a bright blue companion 
may be detectable. 

In order to constrain the contribution of a BB compo¬ 
nent in our 2014 spectrum of SN 2011fe, we need a template 
-1-981 day nebular spectrum that is dominated by SN la flux 
to which we can add synthesised BB spectra. Since we do 
not yet have such a late-time spectrum for any other normal 
SNIa, we instead use the 2013 spectrum, flux scaled to the 
2014 spectrum as described in (We do not use the synthe¬ 
sised 2012-1-2013 composite spectrum for this first method 
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Figure 5. The 2014 spectrum of SN 2011fe (black) and the flux- 
scaled 2013 spectrum (blue; scaling described in § |2.3[ l compared 
with blackbody spectra of temperatures and luminosities in the 
range predicted for a PIRS from Pan et al. (2012, 2013). 


Figure 7. The mean residual, H, in the 4480—5500 A range. Bars 
represent one standard deviation on the distribution of residual 
values (e.g.. Figure]^. Horizontal dotted lines are drawn at 7^ = 
0, 1, 2, and 3. It is evident how the temperature of the PIRS is 
not as constraining as the luminosity. 


to constrain the BB component because it is sensitive to 
the flux calibration and changes in relative features, but we 
do use it below in Section 3.2.11. We then synthesise a set 
of BB spectra, and add them to the 2013 spectrum. To il¬ 
lustrate this idea, the 2013 spectrum is shown in Figure 
along with our 2014 spectrum of SN 201 Ife, and BB spectra 
with the range of temperatures and luminosities predicted 
for PIRS from Pan et al. (2013), to which our test is senstive: 
log(refi/K) = 4.0-4.6 and log(I//L©) = 3.0-4.0. 

We compare this set of synthesised SN-I-BB spectra with 
the “Redux A” spectrum of SN 2011fe from 2014 because 
it shows an excess of blue flux over “Redux B” (!j^, and 
will produce more conservative upper limits on the PIRS 
BB parameters. To statistically assess at what temperature 
and luminosity the SN-I-BB spectrum is significantly dif¬ 
ferent from the 2014 spectrum of SN 2011fe, we first flux- 
match the SN-I-BB spectrum to the 2014 spectrum on the 
red side (6000-7000 A). On the blue side, we then subtract 
the smoothed 2014 spectrum (/sn) from the SN-I-BB spec¬ 
trum (/sn-i-bb) and divide this by the error spectrum of the 
2014 data (bsn) to obtain a residual spectrum, 


^ /sN-rBB(A) - /sn(A) 
esN 

In Figure we show examples of our SN-I-BB models com¬ 
pared with the 2014 spectrum of SN 2011fe, and the resulting 
distribution of residuals, TZ, at the blue end of the spectrum 
(4480-5500 A). When the mean and standard deviation of 
the residuals are inconsistent with zero, we can rule out the 
presence of a PIRS with the luminosity and temperature of 
the BB model added to our 2013 SN spectrum. This method 
is somewhat similar to the work of Fox et al. (2014), who 
search for the signature of a companion star in a spectrum 
ofSN1993J. 

In Figure we plot the mean residual, TZ, for the vari¬ 
ety of PIRS considered. We find that temperature is better 
constrained than luminosity, and that we can rule out the 
presence of a PIRS with T < 10000 K and L > 10000 L© at 
the 2cr level. As seen in Figure within the spectral range 


of 4480-5500 A this test for a large change in the continuum 
is mostly insensitive to the small changes in the features be¬ 
tween epochs - but not entirely. The mean residual does 
not go to zero when the lowest-luminosity BBs are added 
because in the range of 4480-5500 A the 2014 spectrum’s 
flux has declined with respect to the neighbouring broad 
feature at ~ 5400 A, as seen in Figure Thus, our 2014 
spectrum is sensitive only to the coolest and brightest BB 
components that can make an undeniably large change in 
the optical spectral continuum, over and above the existing 
small differences in features. 

Compared to the potential nondegenerate He PIRS 
models of Pan et al. (2013), our temperature and luminosity 
limits are not very constraining on the type of companion 
star at 3 yr. However, all model PIRS luminosity evolution 
tracks peak at ~ 10 yr, and furthermore the full parameter 
space for all types of companions is not covered by these 
models; future theoretical work may find that other types of 
stars may be cooler and brighter. With patience we still have 
a chance, and there is some recent encouragement in this re¬ 
gard. Post-explosion imaging of SN lax 2008ha obtained 4 yr 
after peak brightness shows a source, sufficiently nearby to 
be a PIRS from the progenitor system, that is L ~ 10'* L© 
(Foley et al. 2014). Future attempts to constrain the pres¬ 
ence of a PIRS will have the benefit of our experience, and 
will be able to use this ~ 1000 day spectrum of SN 201 Ife 
as their template. 


3.8.1 Variance Smoothing by a BB Component 

A second way to constrain the presence of a PIRS compan¬ 
ion star in the 2014 spectrum of SN 201 Ife is to compare the 
amplitude of flux variation in 2014 with that of earlier spec¬ 
tra. A significant contribution from a blackbody spectrum 
will “smooth out” the typical features of a nebular-phase 
SNIa. For example, in Figure we compare the smoothed 
spectrum from 2014 with a variety of SN+BB spectra cre¬ 
ated using the 2012-1-2013 composite. These spectra have all 
been renormalised to have an integrated flux equal to that of 
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Figure 6. Left: The original 2014 spectrum of SN2011fe (black, thin) and smoothed (black, thick), and the 2013 SN+BB spectrum 
(blue). Top and bottom rows show the two extrema models of nondegenerate PIRS companions as examples (Tefi and Lq as in plot 
legend). Right: The distributions of residuals {TZ in Equation|^, with the mean (vertical line), standard deviation (horizontal solid line), 
and three standard deviations (horizontal dotted line). 


the 2014 spectrum (i.e., they would produce equivalent pho¬ 
tometric measurements). The amplitude of flux variation is 
clearly suppressed when we add BB components that are 
brighter at optical wavelengths. 

The flux variance, Va„(A), on a velocity scale of An 
km s“^, for a given wavelength A, is expressed by 



EvIa? (/(A')-/a4A))' 


(4) 


where Ai and A 2 are the wavelengths that correspond to 
±At), /At)(A) is the mean spectral flux within ±Au, the nu¬ 
merator is the sum of the residuals between the spectral 
flux and the mean spectral flux within ±Au, and the de¬ 
nominator is the number of pixels within ±Aw. Although 
not exactly the statistical definition of variance, in its way 
V’ai,(A) is a measure of how much the spectrum varies in 
flux in a 2Av km s“^ bin centred on A. As an example, we 
show Vad(A) for Av = 10,000 km s”'^ for the 2012-1-2013 
composite spectrum of SN 201 Ife in Figure 

The value of Vav will also depend on the flux uncer¬ 
tainty, af, or noise of the spectrum. To illustrate this, we 
take the flux error spectrum of the composite 2012-1-2013 
spectrum, multiply it by a constant factor, randomise the 



3500 


4000 4500 5000 

Rest Wavelength [A] 


5500 


Figure 8. To illustrate the effect of adding a BB component to a 
nebular SNIa spectrum, we plot the 2014 spectrum of SN2011fe 
(black, dotted), the 2012+2013 composite spectrum (black, solid), 
and the composite SN+BB for a variety of BB temperatures and 
luminosities (coloured solid lines). All spectra are renormalised 
to have an integrated flux equal to the 2014 observations. The 
temperatures and luminosities of the BB spectra are given in the 
legend at left, and for the reader’s reference we provide velocity 
scale bars in the top-right corner. 
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Figure 9. Flux variance spectra, Va«(A), with a velocity scale 
of Av = 10,000 km s"!, for the SN2011fe composite 2012+2013 
spectrum. To demonstrate the effect of flux-measurement uncer¬ 
tainty (cry) on the value of Va„(A), we add noise to the composite 
2012+2013 spectrum and recalculate the variance. The black line 
represents the original measurement uncertainty; the green and 
red lines represent 2af and 3af, respectively; and the blue line 
represents the addition of flux errors from the 2014 spectrum of 
SN2011fe, cr2oi4. As expected, VadCA) is higher for spectra with 
larger flux-measurement errors. 

sign, add it to the original composite spectrum, and recal¬ 
culate t^Av(A). In Figure]^ the green and red lines show this 
for constant factors of 2 and 3, respectively, while the blue 
line shows this when the flux uncertainty of the 2014 spec¬ 
trum is applied in the same fashion. The 2014 spectrum of 
SN 2011 has the largest flux uncertainty. 

We now use the flux variance to constrain the pres¬ 
ence of a BB component in the 2014 spectrum of SN 2011fe. 
The exact sequence of steps is to (1) compute the BB spec¬ 
trum for a given temperature and luminosity and add it to 
the 2012+2013 composite spectrum for SN 2011fe, (2) renor¬ 
malise such that the SN+BB composite has the same inte¬ 
grated flux as the 2014 spectrum, (3) randomise the sign of 
the flux errors for the 2014 spectrum and add them to the 
renormalised SN+BB composite, and (4) compute the flux 
variance spectrum Vav(A). 

The typical velocity scale of features in SNela is ~ 
10,000 km s“^. In order to show the effect of Av on Va„(A), 
include the results with 1000 and 5000 km s“^ in Figure 
(for the size of Av in wavelength, see the top-right corner 
Figure]^. A value Av = 1000 km s“^ does not sample 
the nebular SNIa features, only the noise, and in Figure 
we can see that the effect of a BB component is not distin¬ 
guishable {Av = 1000 km s“^ is included mainly as a sanity 
check). However, when we use Av = 5000 or 10,000 km s“^, 
we And that the flux variance for the composite SN+BB 
spectrum looks distinctly different from the composite alone. 
While the variance spectra of the 2012+2013 composite and 
the 2014 spectrum are not exactly identical (solid and dot¬ 
ted black lines in Figure |10[ ), qualitatively they are more 
similar to each other than to a SN+BB with T ~ 10000 K 
and log(L/L 0 ) « 4. In this way, we judge that a bright, cool 
FIRS companion is unlikely for SN2011fe. Ultimately, this 
method of flux variance does not provide a tighter constraint 
on a putative FIRS than the flux-residual method described 


in §3.2[ but it is independent of any assumptions regarding 
evolution in the spectral features between the 2012, 2013, 
and 2014 spectra. 


4 CONCLUSIONS 

We demonstrate that extremely late-time photometry and 
spectroscopy of the nearby SNIa 201 Ife shows little evo¬ 
lution between ~ 2 and ~ 3 yr after explosion, and that 
there are several curious aspects of the nebular-phase fea¬ 
tures. There are new features at 5800 and 6400 A that could 
be attributed to [Fel], which may be consistent with the 
material cooling as indicated by the absence of [Felll] and 
[Coin] at these late times. Alternatively, Nal D is a poten¬ 
tial source for the 5800 A feature in a ~ 1000 day spectrum. 
We also consider the possibility that the 6400 A feature is 
the broad, high-velocity hydrogen predicted by Mazzali et 
al. (2014). We hud that a new emission feature has arisen 
between [Fell] A7155 and [Nill] A7370, which could be at¬ 
tributed to [Call] AA7292, 7324. Emission from Na and Ca 
may suggest the presence of circumstellar material from the 
progenitor system. We show that the velocity of the nebular 
features has continued to evolve linearly at ~ 350 km s“^ 
per 100 days, and discuss how the change in line morphol¬ 
ogy suggests that line blending may play a role, which argues 
against the hypothesis that nebular features are exhibiting 
a rest-frame redshift. 

The presence of hydrogen in the CSM would support 
a nondegenerate companion for SN 20 life, but we do not 
see any evidence of narrow (FWHM ~ 1000 km s“^) Ha 
emission in our spectrum. Given this nondetection, we place 
an upper limit on the mass of hydrogen in the system, but 
we find that the models upon which this is based may not 
be appropriate for our extremely late-time spectrum. A lack 
of hydrogen in the progenitor system supports a degenerate 
companion for SN 2011fe. This agrees with the pre-explosion 
archival HST image analysis of SN 2011fe presented by Li et 
al. ( 2011 ), which rules out red giants and most helium-star 
companions. 

In this new, extremely late-time regime of SN la studies, 
the flux from the SN is sufficiently low that a contribution 
from a luminous secondary star could be detected. Although 
all preceding evidence indicates a degenerate companion for 
SN 2011fe, in this paper we make the first observational tests 
for a nondegenerate FIRS. We synthesise a series of SN+BB 
spectra, and statistically constrain the temperature and lu¬ 
minosity for a putative companion with two independent 
techniques based on flux residuals and flux variance. We find 
that only the coolest, brightest FIRS models of Fan et al. 
(2013) are constrained by our observations. Although these 
tests are not currently sensitive enough to provide meaning¬ 
ful constraints on the progenitor scenario, the FIRS model’s 
luminosity evolution does not peak until ~ 10 yr after explo¬ 
sion — so with patience we may yet improve in this regard. 
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Figure 10. Constraining the presence of a BB component in the 
2014 spectrum of SN 2011fe using the flux variance, VavC^)- From 
top to bottom, we show the flux variance when velocity scales of 
Av = 1000, 5000, and 10,000 km s“^ are used. In each plot, 
we show the variance for the 2014 spectrum of SN 2011fe (black 
dotted line), the 2012+2013 composite spectrum with the flux 
uncertainties of 2014 applied (black solid line), and the compos¬ 
ite SN+BB spectra (coloured solid lines, where colour represents 
BB temperature and luminosity as in the plot legend). Variance 
converges at the blue end where smaller-scale features and noise 
dominate. 
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